Casimir forces beyond the proximity approximation 
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The proximity force approximation (PFA) relates the interaction between closely spaced, smoothly 
curved objects to the force between parallel plates. Precision experiments on Casimir forces necessi- 
tate, and spur research on, corrections to the PFA. We use a derivative expansion for gently curved 
surfaces to derive the leading curvature modifications to the PFA. Our methods apply to any ho- 
mogeneous and isotropic materials; here we present results for Dirichlet and Neumann boundary 
conditions and for perfect conductors. A Pade extrapolation constrained by a multipole expansion 
at large distance and our improved expansion at short distances, provides an accurate expression 
for the sphere-plate Casimir force at all separations. 



PACS numbers: 12.20.-m, 44.40.+a, 05.70.Ln 

First introduced by Derjaguin [T], the proximity force 
approximation (PFA) relates forces between gently 
curved objects at close separations to the correspond- 
ing interactions between flat surfaces over an area set 
by the local radii of curvature. Originally developed in 
the context of surface adhesion and colloids, the PFA 
has found application as far afield as the study of inter- 
nuclear forces and fission[21 [S]. In particular, the PFA 
has long been applied to Casimir and other fluctuation 
forces that are short range and therefore dominated by 
the surface proximity. Despite its nearly universal use 
to parameterize Casimir forces in experimentally acces- 
sible geometries like a sphere opposite a plate |U[S], not 
even the flrst correction to the PFA in the ratio of the 
inter-surface separation to the scale of curvature, d/ R, is 
known for realistic dielectric materials or even for generic 
perfect conductors. Casimir forces have non-trivial and 
sometimes counterintuitive dependence on overall shape, 
increasing the interest in developing a thorough under- 
standing of the corrections to the PFA and its range of 
validity. As the experimental precision of Casimir force 
measurements improves, they are becoming sensitive to 
PFA corrections, further motivating renewed interest in 
this subject. For example, an upper bound on the mag- 
nitude of the flrst order correction to the PFA for a gold- 
coated sphere in front of a gold-coated plane was obtained 
a few years ago in an experiment by the Purdue group 
[5], which found that the fractional deviation from the 
PFA in the force gradient is less than 0.4 d/R. 

We here provide a general expression for the flrst cor- 
rection to the PFA for the Casimir energy of objects 
whose surface curvatures are large compared to their 
inter-surface separation. Scattering methods have pro- 
vided a fertile route to computing Casimir forces for non- 



planar shapes [7j. While conceptual threads connected 
to this approach can be traced back to earlier multiple- 
scattering methods fH|, concrete analytical [7] and nu- 
merical [9] results are relatively recent. Indeed, the ex- 
perimentally most relevant set-up of a sphere and plate 
was only treated in 2008 [TD]. The scattering method 
is most powerful at large separation, providing analytic 
expressions for a multipole expansion in the ratio R/d 
and numerical results for R/d ^ 10, but becomes in- 
tractable at the values of d/i? ~ 10~^, needed to study 
the PFA. In a feat of mathematical dexterity, Bordag 
and Nikolaev (BN) summed the scattering series to ob- 
tain the first correction to the PFA for the cylinder/plate 
and sphere/plate geometries, initially for a scalar field 
obeying Dirichlet (D) boundary conditions (bc)|Trj, and 
later for the electromagnetic (EM) field with ideal metal- 
lic bc[Tl]. Intriguingly, whereas the correction to the 
force in the former case for both geometries is analytic in 
d/R, the sphere/plate EM case was predicted to include 
logarithmic corrections (^ {d/R) ln(d/i?)). While the BN 
PFA corrections for cylinders (where there are only an- 
alytic corrections), were independently verified by Teo 
|13) . the results of Ref. [11] [T2] for the experimentally 
interesting sphere/plate EM case remain unconfirmed. 

A recent work by Fosco et. aZ.[n! provides a quite 
promising perspective on Casimir PFA corrections. Rein- 
terpreting perturbative corrections to parallel plate forces 
[TBI [T7] , they propose a gradient expansion in the local 
separation between surfaces for the force between gen- 
tly curved bodies. They implement their program for 
scalar fields and confirm that their general approach re- 
produces the BN cylinder and sphere/plate correction for 
D be. Inspired by this work, the current paper is or- 
ganized as follows: (i) We show that the extension of 
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Ref. [2] to two curved surfaces is constrained by tilt in- 
variance of the reference plane (from which the two sep- 
arations are measured). This provides a stringent test 
of the self-consistency of perturbative results, (ii) Going 
beyond the scalar fields and D be of Ref. [H], we com- 
pute the form of the gradient expansion for Neumann 
(N) , mixed D/N, and EM (perfect metal) boundaries. 
Interestingly, we find that the EM correction must co- 
incide with the sum of D and N corrections, (iii) We 
reproduce previous results for cylinders jl3] with D, N, 
and mixed D/N conditions, and the sphere with D be. 
However, we do not confirm the BN predictions for the 
sphere/plane geometry (12] either with N or EM be. (iv) 
As a further test, we introduce a Fade approximant for 
the (sphere/plate) force at all separations, based on an 
asymptotic expansion at large distances and the gradient 
expansion at proximity. The Fade approximant agrees 
excellently with existing numerical data and implies log- 
arithmic corrections to the FFA beyond the gradient ex- 
pansion performed here (see Eq. (|9|). 

Consider two bodies with gently curved surfaces de- 
scribed by (single- valued) height profiles z = Hi{x) and 
z = H2{x), with respect to a reference plane S, where 
X = (x, y) are cartesian coordinates on S and the z axis is 
normal to S. For our purposes it is not necessary to spec- 
ify the precise nature of the fluctuating quantum field. 
Our results are valid for scalar fields, for the EM field, 
and indeed even for spinor fields. Similarly, the bound- 
ary conditions satisfied by the fields on the surfaces of 
the plates can be either ideal {e.g. D or N for the scalar 
field, or ideal metal in the EM case) or those appropriate 
to a real material with complex dielectric permittivity. 
The only restriction on the be is that they should de- 
scribe homogeneous and isotropic materials, so that the 
energy is invariant under simultaneous translations and 
rotations of the two profiles in the plane E. 

Following Ref. [M], we postulate that the Casimir en- 
ergy, when generalized to two surfaces and to arbitrary 
fields subject to arbitrary be, is a functional E[Hi,H2] of 
the heights Hi and H2 , which has a derivative expansion 

E[Hi,H2]^ [ (fxU{H) [l + l3i{H)VHi-VHi 

JT, 

+ P2{H)VH2 ■ VH2 + /3x {H)VHi ■ VH2 

+ l3^{H)z-iVHiXVH2) + ■■■], (1) 

where H{x) = H2{x)—Hi{x) is the height difference, and 
dots denote higher derivative terms. Here, U{H) is the 
energy per unit area between parallel plates at separation 
H; translation and rotation symmetries in x only permit 
four distinct gradient coefficients, /3i(iJ), (32{H), l3x{H) 
and f3-{H) at lowest order. While the validity of this 
local expansion remains to be verified, it is supported by 
the existence of well established derivative expansions of 
scattering amplitudes [T5] from which the Casimir energy 
can be derived by standard methods [7]. 



Arbitrariness in the choice of S constrains the /3- 
coefficients. Invariance of E under a parallel displace- 
ment of E requires that all the f3's depend only on the 
height difference H, and not on the individual heights 
Hi and H2. The /3-coefficients are further constrained 
by the invariance of E with respect to tilting the refer- 
ence plate E. Under a tilt of E by an infinitesimal angle 
e in the {x, z) plane, the height profiles Hi change by 
SHi = —e(^x + Hi"^^^, and the invariance of E implies 

2{Pi{H) + /32{H)) + 2l3x{H)+ H^^-1 = 0, 
/3_(i/)=0. (2) 

Hence the non-vanishing cross term /3x is determined by 
/3i, /32 and U. 

By exploiting Eqs. ^ we see that, to second order in 
the gradient expansion, the two-surface problem reduces 
to the simpler problem of a single curved surface facing 
a plane, since ?7, /3i and /32 can be determined in that 
case and (3x follows from Eqs. ([2]). Therefore in Eq. ([T]) 
we set Hi = and define l32{H) = (3{H). By a simple 
generalization of Ref. [14 , where the problem is studied 
for a scalar field subjected to D be on both plates, we can 
determine the exact functional dependence of P{H) on H 
by comparing the gradient expansion Eq. ([T]) to a pertur- 
bative expansion of the Casimir energy around flat plates, 
to second order in the deformation. For this purpose, we 
take E to be a planar surface and decompose the height 
of the curved surface as H{x) — d + h{x), where d is 
chosen to be the distance of closest separation. For small 
deformations \h{x)\/d <C 1 we can expand E[0,d + h] as: 

E[0,d+h]= AUid)+fi{d)hiO)+ [ -^G{k;d)\hik)\'' , 

(3) 

where A is the area, k is the in-plane wave-vector and 
h{k) is the Fourier transform of h{x). The kernel G{k; d) 
has been evaluated by several authors, for example in 
Ref.[Tn] for a scalar field fulfilling D or N be on both 
plates, as well for the EM field satisfying ideal metal be 
on both plates. More recently, G{k; d) was evaluated 
in Ref. for the EM field with dielectric be. For a 
deformation with small slope, the Fourier transform h{k) 
is peaked around zero. Since the kernel can be expanded 
at least through order fc^ about k = 0[T11, we define 

G{k;d)^j{d) + S{d)k^ +••• . (4) 

For small h, the coefficients in the derivative expansion 
can be matched with the perturbative result. By ex- 
panding Eq. ([1]) in powers of h(x) and comparing with 
the perturbative expansion to second order in both h{k) 
and fc^, we obtain 

U\d)=fi{d), C/"(d) = 27(d) , /3{d) = ^^, (5) 
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with prime denoting a derivative. Using the above re- 
lation, we computed the coefficient /3 for the following 
five cases: a scalar field obeying D or N be on both 
surfaces, or D be on the curved surface and N be on 
the flat, or vice versa, and for the EM field with ideal 
metal be. Since in all these cases the problem involves 
no other length apart from the separation d, /3 is a pure 
number. /3d was computed in [14) . and found to be 
/3d = 2/3. We find = 2/3(1 - SO/tt^), /3dn = 2/3, 
/3nd = 2/3 — 80/ 77r^ (where the double subscripts denote 
the curved surface and flat surface be respectively), and 
/3em = 2/3 (f - 15/tt^). Upon solving Eqs. ^ one then 
finds /3x = 2 — /3i — /32, where /3i and (32 are chosen to 
be both equal to either /3d, (3n or /3em, for the case of 
identical be on the two surfaces, or rather /3i = /3dn ^nd 
(^2 = /3nd for the case of a scalar field obeying mixed 
ND be. While we obtain agreement for D conditions, our 
results for /3 in the case of N and EM conditions disagree 
with Refs. dH [I2], which finds /3n = 2/3 - S/tt^ and 
/3em = ~2.1, multiplied by logarithmic terms in the lat- 
ter case that we do not obtain. We note that our value 
of /3em is equal to the average of /3d and /3n. Using 



the well known results U{d) 



-an' 



■hc/{lUOd^), where 



ctD — CKN — Q!em/2 = 1 and odn = c^nd = —7/8 in 
Eq. ([T]), it is easy to verify that to second order in the 
gradient expansion the ideal-metal EM Casimir energy 
for two arbitrary surfaces is equal to the sum of the D 
and N Casimir energies in the same geometry. 

Using the values for /3 and /3x obtained above, it is 
possible to evaluate the leading correction to PEA, by 
explicit evaluation of Eq. ([T]) for the desired profiles. For 
example, for two spheres of radii Ri and i?2, both with 
the same be for simplicity, we obtain: 



E — EpFA 



where £^pfa 



d 



+ (2/3-f) — + 



d 

Ri 



- (ttTT^fici?! i?2 ))/ [1440^2 



d 

" (6) 
i?2)]. The 

corresponding formula for the sphere-plate case can be 
obtained by taking one of the two radii to infinity. Sim- 
ilarly, for two parallel circular cylinders with different 
radii and our results, not given here for brevity, fully 
agree with those derived very recently by Teo [T3]. Fi- 
nally we consider two circular cylinders whose axes are 
inclined at a relative angle 0, 



E : 



aiT^hc^ R1R2 
720dsin6' 



1 



8 1 Ri + R2 



(7) 



again assuming the same be on both surfaces. 

Surprisingly, a hyperboloid/plate configuration, 
h{r) = ^/ R'^ + A^r^ — R leads to a PEA correction pro- 
portional to (2/3 -I- 1/X^)d/R, which vanishes for the EM 
case if A = f.20. To better understand what geometric 
features of the surface determine the PEA corrections, 
we considered a general profile, expanded around the 
point of shortest separation, including cubic and quartic 



terms in addition to quadratic terms representing its 
curvature. The cubic terms - representing an asymmetry 
absent in the shapes considered above - actually result 
in a correction term that scales as ^/d/R. The more 
standard correction proportional to d/R depends on 
both cubic and quartic but no higher order terms in the 
expansion of the profile. The quartic coefficients in fact 
have opposite signs for the sphere and hyperboloid, and 
it is this difference that accounts for the possibility of a 
vanishing PEA correction (for the hyperboloid with EM 
be). 




FIG. 1: Casimir energy normalized to PFA energy as a func- 
tion of d/R for the sphere-plate configuration: Numerical data 
(dots, from Ref. [10,), Fade approximants (solid curves) and 
the first correction to the PFA (dashed lines). Also shown 
are the linear slopes (dashed-dotted lines) computed by BN 

[niis]- 

The disagreement between our results for N and EM 
be and those of BN[TT1 [T^] in the prototype configura- 
tion of sphere and plate calls for a comparison with in- 
dependent results. Precise numerical results based on 
an extrapolation of a spherical multipole expansion ex- 
ist for d/R > 0.1 [TU]. Also, at large separations ana- 
lytic asymptotic expansions (AE) of the rescaled force 
/ = R^F/{hc) ~ X^J^i fj''"''"'^'' ill powers of r = R/d are 
available for sufficiently large n, where jo = 2 for D and 
jo = 4 for N and EM conditions, respectively. We resum 
the AE for /, match it with our improved PFA result 
using the method of Pade approximants, and compare 
the resulting interpolation with the numerical results of 
Ref. [To]. We introduce a Fade approximant f[M/N]{r) 
that agrees with n terms of the AE for small r and with 
the leading two terms ~ r^, for large r given by Eq. (|6|. 
This requires /V = Af — 3 whereas M is determined by 
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n, so that all coefRcients in 



flM/M-3] (r) 



Pq + Pir + P2r + . . . + Pmt 



M 



1 



qir + q2r- 



M-3 



(8) 



are uniquely determined. The energy, obtained from 
/[M/M-3] by integration over distance is shown with 
n — 7, 13 and 15 for D, N and EM be, respectively, in 
Fig. [l] together with our result of Eq. ^ for the sphere- 
plate case (i?i = i?, i?2 — oo) and the numerical re- 
sults of Ref . [To] . The difference between Fade approxi- 
mants and the numerical data varies between 0.009% and 
0.062% (D), 0.048% and 0.42% (N), 0.009% and 0.248% 
(EM) . This agreement is quite remarkable given that the 
Fade approximants arc obtained independent of the nu- 
merical data. The Fade approach may thus be valuable 
for Casimir interactions between slowly varying surfaces 
at all separations starting from expansions valid at small 
and large distances respectively. Figure [T] also shows the 
BN results for the first correction to the FFA. 

The Fade approximant for the force at small d/R and 
integration yields for the energy the expansion 



E 



EpFA 



1 



R 



.2,^. l0g-+0 




(9) 



with 6*1 = 2/3 - 1 and 6*2 = 0.08 for D, 02 = -24.01 for 
N and 62 — —4.52 for EM conditions. The appearance 
of a logarithm at second order in Eq. (|9| follows from 
our assumption that the force can be expanded in pow- 
ers of d at small d including a term ^ l/d[20 . A Fade 
approximation to the energy would yield no 1/d term in 
the force. To check that the logd term in the energy 
is necessary, we have tried a Fade analysis of the energy, 
which of course contains no logarithms, and found much 
poorer agreement with the numerical data (with frac- 
tional differences of the order of 100 times larger than 
those given above for the Fade approximant for the force 
and moreover a pole in the approximant). One could 
perform Fade approximations to derivatives of the force 
which might turn out to be even more accurate but this 
would only change the functional form (logarithmic fac- 
tors) of Eq. ([9| at third order in d/R. As an independent 
test of our result 61 = 2/3—1 we fitted the functional form 
of Eq. ([9]) to the numerical data reported in Ref. [ini HH] ■ 
The corresponding values for 0i are summarized in Table 
|lj When comparing numerical and analytic values for 61 
one must consider that they are limited to d/R }^ 0.1 
with somewhat better accuracy for the D and EM case. 
Finally, we note that the presence of the logarithm at or- 
der {d/RY indicates that the gradient expansion for the 
energy might break down at next order. 

Based on the derivative expansion proposed by Fosco 
et al. [15, we have described a systematic treatment of 
the first non-trivial corrections to the FFA. While this 
paper focused on ideal boundary conditions, the scatter- 
ing amplitudes for material with arbitrary EM response 



BOUNDARY 


6*1 FROM FIT TO 


9i FROM GRADIENT 


CONDITION 


NUMERICAL DATA [TOj 


EXPANSION 


D 


-hO.36 


1/3 


N 


-2.99 


l/3-40/7r2 = -3.72 


EM 


-1.62 


l/3-20/7r2 = -1.69 



TABLE I: Coefficient 01 for D, N and EM boundary con- 
ditions on sphere/plate from a fit of the numerical data of 
Ref. |10j to Eq. (|9| and from the gradient expansion. 

(frequency dependent e(w),/i(w), isotropic or not, local 
or not) have been computed in the gradient expansion 
and can be used to compute improved FFA corrections 
in experimentally relevant situations of sphere/plate as 
well as corrugated plates and other setups. 
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